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Abstract. Rapid urbanization and concretization are the main sources of formation and 
existence of Urban Heat Island (UHI). Due to high concentration of pollutants in urban 
environments, the residents are exposed to unexpected health issues. This study aims at 
delineating the temporal variations in the spatial extent of UHI over Gujranwala using Landsat 
thermal imagery. It also aims at determining the variations in pollutant concentration in the 
atmosphere due to vehicle's tailpipe emissions and fossil fuel burning by industrial plants. We 
used various indices e.g., NDVI, NDBI, NDWI and land surface temperature calculations to 
investigate spatiotemporal variations in urban growth patterns and their impacts on the UHI. 
The results show that the UHI enlarged in all direction specially in the north west during the 
study period which is similar to urban growth trends. There also exist positive correlation 
between industrial and vehicle's discharge with pollutant’s concentration in atmosphere. 
Remote sensing tools proved elegant in trend mapping and analysis. 


Keywords: Urban heat island, Normalized Difference Vegetation Index (NDVI), Normalized 
Difference Built-up Index (NDBI), Normalized Difference Water Index (NDWI), Land 
Surface Temperature (LST). 
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1. Introduction. 


The recent migration trends of 
masses toward developed cities in search of 
basic needs of life is very common around 
the world. It is the main reason of extensive 
urbanization [1] which has transformed the 
metropolitan areas into heat islands [2]and 
this phenomenon is termed as Urban Heat 
Island (UHI) [3]. The speedy urbanization 
has enhanced the anthropogenic activities 
that have led to the development of urban 
microclimates. Microclimates are getting 
warmer due to many factors including 
tailpipe emissions by vehicles used for 
transportation, fossil fuel and coal burning 
by industrial plants and the cooking 
concentration at domestic and commercial 
level [1] The growth of vehicles adds 
greenhouse gasses and (Chloro Fluoro 
Carbons (CFCs) in the atmosphere which 
absotb the outgoing long wave radiations 
that contribute to increase the city 
temperature more than periphery. Urban 
heat islands are of two types: 1) Atmosphere 
heat island 2) Surface heat island. The 
atmospheric heat island is further subdivided 
in two types: 1) Urban canopy layer; and 2) 
Urban boundary layer. 

The surface heat island can be 
studied in detail due to availability of freely 
available remotely sensed datasets [3]. These 
datasets have motivated researchers to 
investigate the main sources of urban heat 
reservoirs currently available worldwide 
[4].The urban heat can be analyzed by 
extracting the Land Surface Temperature 
(LST) through satellite derived data. A wide 
range of thermal datasets recorded by various 
sensors including Advanced Spaceborne 
Thermal Emission and reflection (ASTER), 
Advanced Very High Radiometer (AVHRR), 
Moderate Resolution Imaging Spectrometer 
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(MODIS), Landsat 5 Thematic Mapper 
(IM),  Landsat-7/Enhanced Thematic 
Mapper (E TM) and Landsat-8 (TIRS 1,2) are 
useful to extract the pixel-based surface 
temperature. Various researches have 
focused on the relationship of Land Use 
Land Cover (LULC) with LST [5],and found 
that the changes in LST are dependent upon 
the spatial arrangement of urban settlements. 


The natural vegetation is the main 
source of carbon sinks which is removed to 
established new settlements [6, 7]. The 
negative impacts of new settlements are 
extensively studied to determine the 
influence of LULC on urban environment 
e.g. Ding et al, 2013 [8] investigated the 
relation of LULC with LST using Landsat 
images. Similarly, Jusuf et la., 2007 [9] used 
LST data to explore the consequence of 
LULC changes on the spatiotemporal 
pattern of UHI in Singapore. Zang et al. 2013 
[10] evaluated the rate of increasing 
population on UHI patterns in Shanghai. Li 
et al, 2011 [11] computed LST and 
Normalized Difference Vegetation Index 
(NDVI) to understand the impacts of urban 
sprawl on the landscape of Shanghai. Kuang 
et al., 2015 [12] examined the difference of 
LST for different LULC classes and pointed 
out that urban areas were found warmer than 
vegetative areas. The thermal properties of 
land have been altered by modification in 
LUCL of an area that causes the change in 
energy budget, creating the UHI. Zhang 
2008 [13] studied various urban factors 
including urban area, water body, vegetative 
area, and the mean NDVI to map the spatial 
extent of UHI. Singh 2004 [14] monitored 
the seasonal variations in LST of Dehli using 
Landsat TM data. NDVI is a common 
indicator used to discriminate the vegetation 
existing in an area from other features like 
built-up area and the water body [15]. Lo et 
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al., 2003 [16] used NOAA-AVHRR data to 
compute temporal variations in NDVI and 
LST to map LULC changes in Iberia. Similar 
studies were conducted by Kawashina 1994 
[17] for Tokyo, the pearl river delta China 
[18] , the Atlanta metropolitan area [15, 19, 
20] through Landsat thermal datasets. Many 
other studies explored the relationship 
between LST, land surface emissivity and 
NDVI [21, 22, 23, 24] 


The main objectives of this research 
were to delineate the temporal variations in 
the spatial extent of urban heat island over 
Gujranwala and to appraise the impacts of 
various pollutants involved in heat traps 
using satellite images from (1995-2016) in 
ARC GIS interface. 

2. Materials and methods. 
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2.1 Investigation site. 


Gujranwala is among the seven most 
populous areas in Pakistan. It is located in 
northeast of Punjab province in Pakistan. It 
is bordered by the districts including 
Sheikhupura in the south, Hafizabad, and 
Gujrat in west and Sialkot in the northeast. 
Gujranwala is known as a third largest 
industrial center after Karachi and Faisalabad 
in Pakistan that contributes to 5% in national 
Gross Domestic Product (GDP) (Gujranwala 
Wikipedia 
https://en.wikipedia.org/wiki/Gujranwala). 
It has hot and semi-arid climate where 
temperature ranges between 36 to 46 °C in 
summer and drops to 5 °C in winter. The 
spatial extent of Gujranwala is mapped in 
Figure 1. 
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Figure 1. Study site. 


2.2 Material and methods. 


February 2019 | Vol 1 |Issue 1 


The spatiotemporal patterns of LST 
are explored and compared using various 
indices including NDVI, Normalized 
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Difference Built-up Index (NDBI) and the 
Normalized Difference Water Index 
(NDWI). The contemporary advantage of 
this research was to use remote sensing data 
in wider issues relating to the impacts of 
variations in NDVI and LST on the 
environment and human health. We used 
Landsat satellite imagery listed in Table 1, 
freely available on United State Geological 
Survey (USGS) website 
(https://earthexplorer.usgs.gov) to estimate 
the spatiotemporal variations in vegetation, 
built up and water indices. Thermal datasets 
are the valuable products of Landsat series 
that return pixel-based temperature values. 


Table 1. Landsat image acquisition dates. 


St No Date Landsat Series 
1 June18, 1995 Landsat 5 
2 May 06, 2000 Landsat 7 
3 June 21, 2005 Landsat 5 
4 June 11, 2010 Landsat 5 
5 May 26, 2016 Landsat 8 
2.3 Pixel-based temperature 
calculation for Landsat 8 thermal 


bands. 


Landsat 8 started functioning on 
April 10, 2013, with a wider range of spectral 
bands within visible, infrared and thermal 
wavelengths. Band 10 (B10) and Band 11 
(B11) are known as thermal bands commonly 
used to compute pixel-based temperature 
values. A thermal dataset is comprised of an 
array of pixels and each pixel preserves a 
unique number called Digital No (DN). 
Initially, the DN values are converted to 
irradiance using Equation 1 as below, 


Irradiance = (3.342x10* x Thermal 
band)--0.1 (1) 
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The value 3.342x10^ is mentioned in 
the metadata file (.M'TL) saved by a complete 
Landsat image. Equation 1 is useful to extract 
pixel-based irradiance in (watt/m'srad). The 
irradiance-based dataset is further converted 
to temperature using the Equation 2 as 
below, 

K 


—— ~) 213 
In(éK, / Irradiance+ 1) 


(2) 
Where Kı and K; ate constants and their 


values are mentioned in the metadata file of 
each Landsat 8 image as below. 


Tet 


Kı for B10 = 774.89 
Kı for B11 = 480.89 
K> for B10 = 1321.08 
Kfor B11 = 1201.44 


Where e is emissivity and its value is taken as 
0.95 [25]. The complete procedure of 
extraction of pixel-based temperature values 
is mentioned on USGS website. 
(http:/ /landsat.usgs. 

gov/Landsat8 Using Product.php). 


2.4 Pixel based 
calculation for 
thermal bands. 


temperature 
Landsat 5,7 


Thematic Mapper (TM) in Landsat 5 
and Enhanced Thematic Mapper Plus 
(ETM+) in Landsat 7 also offer thermal band 
B6. However, the procedure to analyze pixel- 
based temperature variations is quite similar 
for the both thermal bands produced by TM 
and ETM- [26]. Initially, the DN values are 
converted to spectral radiance using the 
Equation 3 as below, 


L- Lait (Las Lais) X DN/255 
(3) 
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Where L, is spectral radiance, Lminand Lmaxate 
the variations in radiance. The values 
assigned to Lminand Lmaxate 607.7 and 
1260.56 respectively. We used the Equation 
4 to extract pixel based temperature values in 
centigrade 


T= Ga 2273 
In(K, / L, 4-1) 
(4) 


2.5 Calculation of built up, vegetation 
and water indices. 


NDVI is a vegetation index 
commonly used to determine the spatial 
extent and health of vegetation in an area. It 
is widely used index that considers a spectral 
band within near infrared wavelength and the 
other in red wavelength using the Equation 5 
below [4, 18, 20]. 

NDVI = NIR -Red (5) 
NIR + Red 

NDVI ranges between -1 to +1 
where +1 highlights the vegetation and -1 
represents all the feature other than 
vegetation. 


NDBI is commonly used index to 
map the built-up areas and barren lands. 
Built-up areas experience a drastic increase in 
reflectance in the Shortwave Infrared 
(SWIR) wavelength range (1.55-1.75)um and 
are markedly lower in the Near Infrared 
(NIR) wavelength i.e., (0.76-0.90)um. The 
differentiation between these two 
wavelength ranges results in positive values 
for built-up pixels [27] as follows, 


| SWIR — NIR 
SWIR + NIR 


NDWI is the water index that is 
commonly used to enhance the reflectance of 


NDBI (6) 
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watet bodies within a study site. The formula 
used to compute NDWI 1s as follows, 


NDWI = Green — NIR 7) 
Green+ NIR 


Water reflects maximum in green 
wavelength range (0.52-0.60)um and absorbs 
all the heat radiations within NIR range. 


2.6 Regression analysis. 


Regression is a statistical indicator 
used to determine the strength of the 
relationship between two variables. We used 
regression to compute the relationship for 
spatiotemporal variations between 
concentrations of various pollutants in 
comparison to their sources. 


2.7 GIOVANNI. 


GIOVAANI is web-based forum 
sponsored by National Aeronautics and 
Space Administration (NASA) that provides 
spatiotemporal variations in concentration of 
atmospheric gases globally. We used 
GIOVANNI to extract the information 
about atmospheric constituents that play a 
vital role in increasing the LST by capturing 
the terrestrial heat. 


Results and discussion. 


We selected the satellite imagery for 
the months May and June to compute the 
vegetation indices e.g., NDVI, NDWI, and 
NDBI for the study site. During these 
months the study area receives maximum 
sunshine hours and peak solar flux which is 
important to delineate UHI. Also, we can get 
high contrast for the indices to map the 
extent of major features including built-up 
area, water body, and the vegetation. NDVI 
is an important indicator commonly used to 
map the spatial extent of vegetation existing 
in an area. We used Equation (5) to generate 
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NDVI based datasets for all satellite images 
mentioned in Table 1 using spectral bands in 
NIR and red wavelength ranges and mapped 
the results in Figure (2) A1-A5, where the 
lush green areas are vegetation and the red 
regions are non-vegetative features. NDBI is 
an elegant indicator that highlights the built- 
up features including buildings, road 
networks and open areas etc. This index is 
used to distinguish urban areas in 
comparison to other features like vegetation 
or watet body. We used Equation (6) to 
generate NDBI based datasets for all Landsat 
images mentioned in Table (1) using the 
spectral bands in SWIR and NIR wavelength 
ranges and mapped the results in Figure (2) 
(B1-B5). NDWI is an index which is used to 
make prominent the water body existing in a 
dataset. Water is a good absorber of heat in 
infrared wavelength that reflects nothing 
toward satellite's sensor, therefore, we get 
waterbody in blackish tone in a remotely 
sensed dataset hence water area is difficult to 
demarcate properly. NDWI is used as a 
substitute to delineate water features existing 
in an area. We applied Equation (7) to all the 
datasets mentioned in Table 1 to delineate 
the water body and mapped results in Figure 
2 (C1-C5). 


Temperature effects the human 
activities and consumption of biomass of an 
area. Various crop promoting activities to 
achieve targeted productivity level is 
operated by temperature. These activities 
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include 1) the impact assessment of net 
radiations on plant growth 2) the effects of 
elevated temperature on biomass generation 
and 3) the effects of heat stress on the water 
intake capacity by plants, therefore, it is 
important to manage the temperature range 
which is considered optimum for the proper 
growth of a particular plant. 


The heat generation in an area is 
largely dependent upon various factors 
including the existing land use, number of 
vehicles, industrial emissions, and burning of 
fuels. The temperature of the study site is 
increasing every year that needs to be 
monitored to create a balance between the 
utban area in comparison to its surroundings. 
We computed pixel-based temperature 
datasets using thermal bands for all satellite 
images mentioned in Table 1 and mapped the 
results in Figure 2 (E1 to E5). We extracted 
vegetation from all NDVI based datasets, 
built-up areas from NDBI datasets, water 
from NDWI based datasets and made 
composites by taking their integrated impact 
e.g., the composite made for the year 1995 
consisted of the collective impact of three 
indices using raster calculator utility in Arc 
GIS 10.1. This method of generation of 
composites is comparatively accurate as per 
other methods of feature mapping using 
classification utility. We used Equations (1,2) 
for thermal datasets of Landsat 8 and the 
Equations (3,4) for Landsat 5,7 to compute 
pixel-based temperature values. 
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Figure 2. Spatiotemporal variations in NDVI, NDBI, NDWI and LST (°C). 
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Spatiotemporal variations mapped in Figure 
2 determine that the total area under 
investigation was 861 km’. This area was 
partitioned by three dominating features 
including built-up area, water body, and the 
vegetation. The maps (A1 to A5) in Figure 2 
are showing a decline to the vegetative area 
of 725 km^was under vegetation in the year 
1995 that degraded to320 km? in 2016. 
Whereas the maps (B1 to B5) in Figure 2 
desctibes progressive jumps to built-up area. 
Built up area was estimated as 93 km? in 1995 
that increased up to 318 km^in 2016. Built- 
up area was observed with a remarkable 
boost of 78 km? during 2010 to 2016. The 
maps (C1 to C5) in Figure2 highlights the 
water body which seems with no change in 
area under water flow. 'The water area was 
observed as43 km’ in 1995 that remained 
unchanged in 2016. We observed a slight 
increase in water area in the map C4 in Figure 
2 which was the accumulation of rainfall 
water that removed after some days by 
evaporation and seepage by the earth. The 
composites containing an amalgam of all 
existing features were mapped (D1 to D5) in 
Figure 2 which show a degradation to 
vegetation in comparison to built-up areas. 
Spatiotemporal variations in land use 
patterns (D1 to D5) were analyzed that 
determine the temporal shift of built-up areas 
toward north-west direction at large scale 
and a reasonable sprawl in all directions at 
small scale. Maps (E1 to E5) in Figure 2 
determines the spatial variability of lands 
surface temperature. The comparison of 
maps (D1 to D5) with (E1 to E5) in Figure 2 
shows the same trend for hotspots marked 
with red color which means that the spatial 
extent of heat islands is highly dependent on 
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built-up land use patterns, e.g., the shape of 
heat island in map E1 was centralized 
according to the built-up extent mapped in 
D1 that show similar patterns for heat island 
as pet temporal changes in built-up land use. 


2.8 Effects of atmospheric gasses on 
urban heat islands. 


The concentration of various gases 
e.g., (CO», CO, NO, and SO;) plays a vital 
role in increasing the atmospheric 
temperature. These gases are capable enough 
to capture longwave radiations emitted by 
the surface of earth creating a natural 
greenhouse effect. The impacts of these 
gasses cannot be ignored while investigating 
the UHI phenomenon in an industrial region 
like Gujranwala therefore, we collected 
temporal variability of these gasses from 
GIOVANNI for our study area and analyzed 
the main sources of production of air 
pollutants in term of CO», CO, NO»and SO». 
These pollutants are highly involved in 
blockage of urban heat and the terrestrial 
radiations which cause to make the urban 
ateas as heat islands. There are many sources 
which produce these pollutants but there are 
two major sources including the fuel burning 
by the vehicles and the industrial emissions. 
These emissions cause to increase the 
concentration of these gases in atmosphere 
therefore, we collected all the data regarding 
the temporal increase in a number of 
industries and the vehicles in Gujranwala 
from Pakistan Bureau of Statistics annual 
reports for the years 1995 to 2016 as 
mentioned in Table 2. 
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Table 2. Temporal variations in concentrations of CO», CO, NO» and SO» from 1995 to 2016. 


Year CO CO; NO; SO; Vehicles Industries 
ppbv PPM 1cm” (x10") — Kgm? (x10°) 

1995 105 363 4.2 2.6 550 
1996 105 365 4.5 2.9 611 
1997 104 368 4.7 3.4 684 
1998 115 370 4.8 3.1 710 
1999 114 371 5 4.6 695 
2000 112 373 5.2 4 25730 695 
2001 105 375 5.5 4.1 122659 620 
2002 110 376 5 4.13 132224 641 
2003 110 378 5.4 4.15 141789 643 
2004 111 379 5.5 4.2 158795 646 
2005 109 380 5.6 4.3 183877 646 
2006 113 382 5.7 4.4 213771 714 
2007 114 383 5.8 4.7 247029 817 
2008 114 386 6 4.5 283154 949 
2009 115 389 6.3 5 327533 1039 
2010 115 392 7 4.6 372733 1056 
2011 116 394 8.1 5 422902 1079 
2012 117 396 8.6 5.2 474825 1092 
2013 118 397 8.8 5.4 534905 1105 
2014 121 404 9.1 5.8 640423 1218 
2015 124 403 9.3 5.7 689294 1278 
2016 130 420 11.5 6.5 796751 2131 


We generated a linear regression to determine the strength of relationship of industrial 
emissions with atmospheric pollutants and mapped the results in Figure 3. 
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Figure 3. Linear regression model applied to various gasses (CO2, CO, NO» and SO) effecting 


the urban temperatures. 


Carbon monoxide (CO) is a 
colorless, odorless and poisonous pollutant 
which is formed by fossil fuel combustions 
emitted by the both vehicles and the 
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industrial emissions. The number of 
registered industries in Gujranwala was 
increased up to 2131 in the year 2016 which 
was only 550 in 1995. CO concentration in 
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the atmosphere was largely influenced by 
industrial emissions as there were 740 
industries in 1998 and the CO concentration 
was 115ppbv as recorded by GEOVAANI 
that reduced to 110ppbv against 641 
industries in 2002. We applied a linear 
regression model to the increasing 
concentration of CO in comparison with 
temporal variation in a number of industties 
and the vehicles to determine the strength of 
relationships. We found that CO 
concentration in the atmosphere was largely 
influenced by the number of vehicles with 
R’*=0.919 as compared to industries as 
mentioned in Figure 3 (a,b). CO blocks the 
heat and creates a natural heat island. 


Carbon Dioxide (CO2) is emitted by 
tailpipe from trucks, cars, vehicles, and 
buses. As the number of vehicles increases 
the concentrations of CO» increases in the 
atmosphere which causes the main source of 
heat trap in the atmosphere. The burning of 
fuels by vehicles add a high concentration of 
CO»1n the atmosphere as determined by the 
regression analysis applied in Figure 3 (d) 
with R? = 0.96. However, industries added 
healthy quantity of CO» in the atmosphere 
analyzed using regression analysis of CO» 
with a number of industries that generated 
R/-0.85 in Figure 3 (c). This large 
concentration of CO» in the atmosphere is 
capable to block terrestrial heat and to create 
a heat island. 


Nitrogen Dioxide (NO; is a 
pollutant which plays a vital role in in making 
ground-level ozone and particulate matter. It 
causes the human lungs infection and 
destroys the immunity system which works 
in defensive mode against respiratory 


infections such as influenza and Pneumonia. 
Sulphur Dioxide (SO;) is a very harmful 
pollutant for the human body that is 
produced by the burning of diesel and coal. 
The existence of SO? in atmosphere is the 
clear indication of working of industrial 
plants. It is highly risky to intake by young 
ones that cause asthmatic attacks. We 
generated a correlation between temporal 
variations of SO»; concentration in the 
atmosphere in comparison to the increasing 
number of vehicles and industries and found 
a strong relationship for the both as 
mentioned in Figure 3 (g,h). 


3. Conclusion. 


Anthropogenic activities played a 
vital role in enhancing the concentration of 
atmospheric pollutants in the atmosphere of 
study site which has created an alarming 
situation not only for human survival but 
also for agricultural productivity. The 
atmosphere of Gujranwala is getting heated 
and this overheating has influenced the crop 
phenology and a seasonal shift in local 
climatic events, therefore, we get rainfall in 
Gujranwala at the time which is not suitable 
for crops e.g., in ripening periods. On the 
other hand, a large number of patients were 
observed in recent decade in outdoors and 
emergencies of hospitals who are hit sun 
stock, dehydration and skin problems. 
People face a number of health issues which 
are due to inhale of these pollutants. CO is a 
precursor of bad ozone and ozone cause skin 
cancer. Therefore, it is important to take bold 
steps at the governmental level to control the 
oxides of Carbon, Sulphur and Nitrogen 
emissions by controlling of fossil fuels and 
coal burning 
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